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We analyze the one-loop vacuum stability and perturbativity bounds on a singlet 
extension of the Standard Model (SM) scalar sector containing a scalar dark matter 
candidate. We show that the presence of the singlet-doublet quartic interaction re- 
laxes the vacuum stability lower bound on the SM Higgs mass as a function of the 
cutoff and lowers the corresponding upper bound based on perturbativity considera- 
tions. We also find that vacuum stability requirements may place a lower bound on 
the singlet dark matter mass for given singlet quartic self coupling, leading to restric- 
tions on the parameter space consistent with the observed relic density. We argue 
that discovery of a light singlet scalar dark matter particle could provide indirect 
information on the singlet quartic self-coupling. 



INTRODUCTION 



The search for the mechanism of electroweak symmetry-breaking (EWSB) is one of the 
primary aims of on-going experiments at the Tevatron and future experiments at the Large 
Hadron CoUider. Until a Standard Model (SM) Higgs boson is discovered at either collider, 
it is by no means assured that EWSB occurs through a low-energy scalar sector that contains 
only a single, complex SU(2)2, doublet. The simplest variant of the latter scenario involves 
the presence of an additional scalar degree of freedom (real or complex) that carries no 
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SM gauge charges. If it exists, a scalar singlet could be the low-energy remnant of a more 
complex higher energy model. Depending on the nature of the scalar potential, a real or 
complex singlet could provide a viable candidate for the observed relic abundance of cold 
dark matter [l lO] and/or give rise to a strong first-order phase transition during the era 
of EWSB as needed for electroweak baryogenesis 11 Its existence could be verified 
experimentally through the discovery of one or more mixed SU(2)l doublet-singlet states 
in colliders |9|, llO| , the observation of a distinctive azimuthal correlation between outgoing 
jets in weak boson fusion | l4-ll6l| . and a signal in the next generation of dark matter direct 
detection experiments jo], IitI. [i8|. 

If a scalar singlet is discovered, a combination of the foregoing experiments could also 
provide information about the parameters of the scalar potential, particularly those that in- 
volve both the doublet and singlet degrees of freedom. Probing the singlet self-interactions 
will be more difficult. This problem is particularly interesting if the scalar singlet provides a 
dark matter candidate, as the nature of the dark sector and its interactions is presently un- 
known. Bounds on the scalar singlet quartic self-coupling can be derived from the dynamics 
of the Bullet cluster [loj and from observations of galactic dark matter halos [l, 0, IM, ISl, 20 



if the scalar is very light {Ms on the order of 1 to 100 MeV) but for heavier masses there 
presently exist no limits. 

It is also interesting to ask, if the SM scalar sector is accompanied by a singlet, at what 
scale A should one expect the high energy desert to be populated with additional degrees 
of freedom? At one-loop order, the absence of a Landau pole in the scalar self-interactions 
should provide a rough guide as to this scale. Since the value of Mh fixes the value of the 
Higgs quartic coupling at an input scale for its renormalization group (RG) evolution, the 
discovery of a heavy Higgs could provide one indication of a maximal value for A. On the 
other hand, if a relatively light scalar is discovered, A cannot in general be too large if the 
vacuum associated with EWSB is to be stable. The corresponding considerations have been 
analyzed extensively in the case of the SM (more recently in Refs. (21I . 22|). with an allowed 



range of 55 GeV ^ Mh ^ (500 — 800) GeV for A ~ 1 TeV and a correspondingly narrow 
range of 150 GeV ^ Mh ^ 200 GeV for A at the GUT scale (neglecting theoretical and 
parametric uncertainties)^. 

To our knowledge, the corresponding expectations for singlet extensions of the SM have 



not been thoroughly studied, although partial analyses were performed in Refs. [23|, |24 
These works considered vacuum stability and perturbativity constraints on a singlet exten- 
sion of the SM in the context of infiation models which identified the infiaton with either 
the Higgs or the scalar singlet. The authors included non-minimal gravitational coupling 
between the Ricci scalar and the infiaton, i.e. the Higgs or singlet, and showed that cosmo- 
logical observables place restrictions on the parameter space which can be stronger than the 
vacuum stability and perturbativity limits. In this work, we do not assume that either the 
Higgs or the singlet serves as the infiaton and so we do not encounter constraints imposed 
by infiation. 

In what follows, we address these issues for the scalar singlet extension of the SM by 
analyzing the one-loop vacuum stability and perturbativity constraints on the theory. For 
simplicity, we restrict our attention to the case of a real scalar singlet S for a potential 



^ As discussed below, the two-loop analysis of Ref. [21] ties the upper bound to the fixed point of the 
two-loop /9-function; for large A, the results are not too different from those obtained using the one-loop 
Landau pole. 
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possessing a Z2 symmetry - henceforth denoted the Z2XSM - and require that the singlet 
vacuum expectation value (vev) (5*) = in order to provide for a stable singlet that may 
be a dark matter candidate. We first find that, in contrast to the SM case, a perturbativity 
criterion based on the avoidance of a Landau pole at large A can lead to an unreliable 
analysis of the vacuum structure. Consequently, a more stringent one-loop perturbativity 
requirement must be imposed. We then show that for a given value of A, the presence of 
the quartic S-H coupling, henceforth denoted 02, lowers both the upper and lower bounds 
on Mh. For given values of and M5, this coupling must generally be sufficiently large in 
order that the relic density of singlet scalars does not over-saturate the observed cold dark 
matter relic density. In this scenario, then, the discovery of a relatively heavy Higgs would 
imply a lower scale for the onset of new physics than in the SM. These considerations are 
particularly pronounced for larger A. 

We also find that requiring that (S*) = be the true vacuum can have interesting impli- 
cations for the singlet scalar self-coupling, denoted 64 below, particularly when consistency 
with the observed relic density is imposed. For given values of M/^, M5, and A, these con- 
siderations can imply a lower bound on the magnitude of this self coupling. Alternately, 
fixing the magnitude of the latter can lead to a lower bound on the singlet dark matter 
mass. To illustrate, for A = 1 TeV and = 120 GeV, a vanishing self-coupling implies 
Ms ^ 50 GeV. Increasing 64 reduces this bound while increasing A raises it. For A > 10^ 
GeV and Mh = 120 GeV, there exists no region of the {Ms, 64, CI2) parameter space that 
can accommodate the relic density. In short, assumptions about the magnitude of the scalar 
self coupling and the scale of new physics can impact the degree to which a scalar singlet 
can explain the DM relic density. 

In the remainder of the paper, we present the analysis leading to these conclusions. In 
Section [21 we review the model potential. Section [3] gives our analysis of the one- loop, RGE- 
improved effective potential and the criteria we choose to maintain a perturbative theory. 
We show illustrative applications of our analysis to the model parameter space and discuss 
the implications of these results in Section |H We summarize our work in Section [51 



SCALAR POTENTIAL 



The most gene ral renormalizable potential involving the SM Higgs doublet H and the 
singlet S is 



6 2 3 4 

where < and where we have eliminated a possible linear term in S" by a constant shift, 
absorbing the resulting S'-independent term in the vacuum energy Q. To set our notation^ 
we take 



V2 



01 + i(P2 

h + Z04 



(2) 



We require that the minimum of the potential occur aX h = 246 GeV = Vq. Fluctuations 
around this vacuum expectation value (vev) are the SM Higgs boson and the 0j are the 



^ For a comparison of our notation and normalizations with those used in the hterature, see the Appendix. 
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would-be Goldstone bosons that become the longitudinal components of the electroweak 
gauge bosons. 

For the case of interest here for which S is stable and may be a dark matter candidate, 
we impose a Z2 symmetry on the model, thereby eliminating the ai and 63 terms. We also 
require that the true vacuum of the theory satisfy (5*) = 0, thereby precluding mixing of S 
and the SM Higgs boson and the existence of cosmologically problematic domain walls. In 
this case, the tree-level minimization conditions give 



^^0 



with masses for the doublet and singlet scalars 



62 + 02 "^0 



(3) 



(4) 
(5) 



to eliminate in favor of Vq. 



where we have utilized the tree-level relation (|3]) 

For purposes of analyzing vacuum stability, inclusion of radiative corrections is essen- 
tial. We begin with the one- loop effective potential, Ves{h, S), whose computation requires 
knowledge of the field- dependent masses for the scalar, fermionic, and gauge boson degrees 
of freedom, denoted mj{h, S). The one-loop contribution, renormalized in MS, is given by 



Vih,S) 



(Hoop) 



E 



647r2 ^ 



ih,S) 



log 



m2 {h, S) 



(6) 



where j runs over all degrees of freedom in the model and where h and S now denote 
the classical fields. For the fermions, we retain only the top quark contribution since the 
contributions from the remaining fermions are suppressed by their Yukawa couplings and 
have a negligible impact on our analysis. For the neutral scalars, allowing for both classical 
fields h and S to be non-vanishing away from the minimum of the potential leads to mixing 
between these degrees of freedom associated with the off-diagonal element in the mass- 
squared matrix 



M\h,S) 



+ f + 02^2 2a2hS 



,2 I Xl,2 

2a2hS 62 + 364^^ + a2h^ 
The corresponding field-dependent eigenvalues are 











f ml 

















{ml + ml) ± ^(^^^^^1)^4^ 



(7) 



(8) 



The remaining field-dependent masses, constant factors Cj, and degeneracy factors rij are 
listed in the Appendix. 

As seen in Eq. (|6]), the one-loop contribution to the potential depends explicitly on the 
arbitrary unphysical renormalization parameter (or 't Hooft scale) /i. In order to ensure 
that the full potential is scale independent, we introduce the running parameters A — )■ 



5 



A {^) , m? — )■ m? (/i) , h ^ h (/i) , etc. that also depend on the scale /i. The full potential 
through one-loop is given by 

Ves{h, S) = V{h, SY''^''^ + V{h, ^)(ii°°P) (9) 

with 

n/^,5)(*-^) = ^/.2 + ^/.^ + ^5V + M/^5^ + ^5^ + fi(/.) , (10) 

where we have suppressed the //-dependence of the fields h = h{ji) and S = S{fi) for 
simplicity . 

The running parameters are obtained from the beta and gamma functions and the cor- 
responding renormalization group equations (RGEs): 

(U) 

= _/^^ ^ d log Zy 

~ Y dfi 2 dfx ^ ^ 

for a given coupling or mass parameter X and field Y [Zy is the wave function renormal- 
ization for Y). The beta functions for the gauge couplings, g, g', and g^, and the top quark 
Yukawa coupling, yt, are unaffected to one- loop order by the addition of the singlet; thus 
they are the same as in the Standard Model and can be found in Ref. 2^. Additionally, 



the wave function renormalizations for the Higgs and the singlet do not receive any one-loop 
scalar contributions. As a result, 7/1 remains the same as in the Standard Model (also in 



Ref. j25|), and 75 = 0. 

All that remains is to find the beta functions for the scalar potential parameters. These 
are determined by requiring scale independence of the effective potential in Eq. to one- 
loop order: 

dV,s{h,S) _ dV{h^y^ dVVhSY^_r, 
H - — /i - h /i ^ — u [L6) 

djj, djj Ojji 

Using Eqs. ( ITOj) and (E]), the definitions for Px and 7y above, and the known results for 7/4 
and 75, we use this scale cancellation to obtain the beta functions for A, 02, &4, m?, and 
62 by comparing terms of the same degree in h and S. The running of the vacuum energy 
is found from terms that are independent of the fields. We have independently verified our 
results, given in the Appendix, by explicitly computing the one-loop MS renormalization of 
the quadratic and quartic terms in the potential and using the scale-independence of the 
bare parameters^. 

As boundary conditions for the running gauge couplings, we take 

g (Mz) = 0.65 
g' (Mz) = 0.358 
gs (Mz) = 1.21 . 



Our results for the beta functions also match those in Refs. 



23 



24|. 
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For the top Yukawa coupling boundary condition, we follow the choice of Ref. 22| : 

rrit (Mt) 



yt 



(Mt) = V2 



where the running top mass, rrit (/i), is related to the top quark pole mass, Mt, by 



M,= (l + l^)m,(M,) 



and the factor ^ (Mt) is explained below. 

For the quartic couplings, we choose the input scale /xq = Mz and will vary \{Mz), 
b^^Mz), and a2{Mz) according to various criteria discussed in Section |3]below. The quadratic 
mass-squared parameter m?{^) will be related to the Higgs boson pole mass using the 
minimization conditions on I4fj. At one-loop order, it is convenient to relate these quantities 
at yU = fo, so we will take Vq as the input scale for both running quadratic parameters^. 

We take as the solution for the Higgs field RGE 

h{fi) = = exp (- [ -fh{ij')di2'/fi')h (14) 

V Jmz J 

with boundary condition 

i{Mz) = \ 

so that 

h {Mz) = h . 

The running of the classical Higgs field h (fi) in f|T^ is very weak - ^ (/x) deviates from one 
by less than 10% over the range Mz ^ fJ^ ^ 10^^ GeV (assuming there are no Landau poles 
present in this range). Meanwhile, S receives no scale-dependence at one- loop order since it 
does not couple to gauge fields and since quartic scalar interactions generate no wavefunction 
renor malizat ion . 

We note that our choice of input scale for the quartic and quadratic parameters would not 
be optimal for a two-loop analysis, since the corresponding one-loop matching corrections 
can be unacceptably large. As discussed in Ref. [2l|, a suitable choice is fiQ = max{Mt, M^}. 
Since we do not encounter any loop level matching conditions in our one-loop Kff analysis, 
the only impact of varying our input scale between Mz and Mt is via the running param- 
eters and the explicit one- loop term in Ves, and the corresponding numerical impact is not 
significant. Consequently, we choose an input scale for each parameter that renders the 
explicit logarithms in Eq. small and that is most convenient for fixing the corresponding 
boundary condition for the running. 

Finally, we comment on a minor difference between our treatment of the scale-dependence 
and that of Ref. j22[ . The authors of that work identified an optimum scale fi* that minimized 
the scale dependence of the effective potential and performed their minimization of the 



^ Note that the beta functions for the couphngs are independent of the 62 and m? parameters, ahowing 
them to be solved independently without the minimization conditions. Once we have determined the 
running of the couplings, and because depends linearly on m'^{fi), we can choose any scale for our 
boundary conditions - vq is the most convenient choice because of the minimization conditions. 
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potential at that scale /i*. In the present work, no such choice of a special scale is made; 
since we include both terms in Qj, the effective potential is scale- independent to one-loop 
order with the residual scale dependence being higher order. When analyzing the behavior of 
Veff at large values of the field as needed for vacuum stabihty considerations, we will choose 
fi so as to minimize the explicit one-loop logarithms, resumming them implicitly through 
the running parameters in the first term of Eq. ([9]) . 



3. VACUUM STABILITY &: PERTURBATIVITY 

It is well-known that in the SM, Ves{h) can develop a minimum at large values of h {i.e., 
h » fo) that is deeper than the electroweak minimum^. This situation arises because of 
the large, negative top Yukawa contribution to /3a that causes A to run negative at large 
scales. In order to avoid the presence of this second, deeper minimum at scales below the 
cutoff A of the effective theory, one requires a sufficiently large value of A at the input scale 
(here taken to be /xq = Mz) so that the turnover of Ves occurs at sufficiently large scales. 
On the other hand, a value of A(/io) that is too large will lead to the presence of a Landau 
pole below A. These considerations lead, respectively, to the vacuum stability and triviality 
bounds on M^, whose value is governed by A at the electroweak scale. 

It is, of course, possible that the SM vacuum is metastable with respect to the second 
minimum at large h with a lifetime longer than the age of the Universe. If during the 
process of electroweak symmetry-breaking the Universe first lands in the finite temperature 
electroweak vacuum, it may remain there as the temperature decreases even if the large 



h minimum eventually becomes deeper. Recent works such as Refs. |32|, |33| have shown 
that allowing for a sufficiently long-lived electroweak vacuum can weaken the lower bounds 
on the Higgs mass by as much as ~ 25 GeV (this change was roughly uniform for cutoff 
scales 10^*^ GeV < A < 10^^ GeV). A careful analysis of the metastability bounds requires 
computing the one-loop corrections to the vacuum tunneling probability, which is beyond 
the scope of the present work; we restrict our attention to the more conservative criterion 
of absolute vacuum stability at T = 0. 

To obtain the one-loop stability bound in the Z2XSM, we first analyze the behavior of 
the potential along the /i-direction {S = 0), and determine the value of X{Mz) such that 
the potential does not develop an absolute minimum for A > h » vq. We choose the 
values of the singlet parameters at the input scale such that no minima occur away from 
the h axis, and subsequently study the impact of departing from this choice. In order to 
obtain the most reliable one-loop analysis, we minimize the impact of large logarithms in 
the explicit one-loop contribution to V^s by choosing = h{Mz) (recall that h{fi) varies 
gently with fi). Since the field-dependent masses along the /i-direction are also proportional 
to h, this procedure leads to the appearance of explicit one-loop logarithms having order 
one arguments, while resumming the large logarithms through the running parameters in 
V{h,S). We then study Ves as a function of h{Mz) = h. An alternate choice, such as 
choosing /i = Mz, would lead to explicit, un-resummed large logarithms for large values of 
h, rendering our one-loop analysis less reliable. 

To obtain the electroweak/dark matter minimum, we choose m?{fi) such that the mini- 



For a comprehensive review, see Ref. [26| and the references therein; also, see Refs. |22l |27H31| 
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mum of the potential occurs for h {vq) ~ h (Mz) = Vq 



dh 



= , (15) 

h=vo, S=0, fJ.=vo 



with positive eigenvalues for the matrix of second derivatives (see below). Except for {vo), 
all quantities in Eq. ( ITSjl are known numerical values (once the RGEs for the couplings and 
the Higgs field have been solved). We subsequently plot Ves as a function of h{Mz), with 
/i = h{Mz) and determine whether or not a second, deeper minimum occurs below A. For 
given values of 02 (M^), 64 (M^), and 62 ('^0); we repeat the procedure for successively smaller 
values of X{Mz) until the appearance of a second, deeper minimum - thereby obtaining the 
minimum value of X{Mz) consistent with absolute stability up to the scale A. 

To arrive at a corresponding upper bound on X{Mz), we begin by fixing hi{Mz), a2{Mz), 
and &2(^^o) and solving for the value of X{Mz) that leads to a Landau pole at the cutoff A. 
We observe that the use of this triviality criterion to establish an upper bound on X{Mz) is 
somewhat subjective. An alternate possibility is to require that the theory remain pertur- 
bative up to the cutoff since the effective potential has been computed using perturbation 
theory. At one-loop order, one minimally must avoid the occurrence of the Landau pole 
below A; doing so, however, does not guarantee that the theory remains perturbative below 
the cutoff. In fact, we find that using this triviality criterion leads to the presence of unsta- 
ble minima at large scales {i.e., large radii in h-S field space) because the couplings in the 
loop corrections become non-perturbative. Whether or not the existence of such minima 
is a physical effect is therefore suspect. To avoid such problems, we choose to adopt more 
stringent perturbative restrictions on the running coupling parameters at the scale A rather 
than rely on the triviality criterion to arrive at an upper bound. 

Detailed studies of the perturbativity requirements in the SM have been carried out using 
next order (two loop) /9-function analyses, resummation techniques, and lattice methods. 
At two- loop order, the Higgs quartic coupling reaches a fixed point, A = App (corresponding 
to Px = 0), rather than a Landau pole. According to the analysis of Ref. [21], the input 
value X{Mz) leading to A(A) = App is significantly smaller than the value that gives a 
Landau pole at A for small values of the cutoff, with the difference decreasing as the cutoff 
is increased. In addition, the work of Ref. [s^] suggests that for A(A) < App/4 the theory 
remains perturbative up to A, while for A(A) < App/2, one may be reaching the limit of 
perturbative behavior. Thus, we expect that for the lower range of our cutoff, use of the one- 
loop triviahty criteria would give the most generous range for the quartic couplings A(M^), 
64 (M^), and 02 (M^:) were it not for the incidence of unstable minima at scales beyond the 
perturbative regime. A full two-loop analysis of the Z2XSM - together with a study of 
scattering amplitudes that depend on these parameters - would lead to a precise statement 
of perturbativity for the quartic couplings and a more restricted range of parameter space. 
However, given the subtle interplay of the various new parameters that occurs already at 
one-loop order and the qualitatively new features that emerge, we defer a refined study of 
two-loop stability and perturbativity considerations to future work. 

In the interim, for the purposes of this analysis we consider two possible criteria to con- 
strain the values of the couplings at the cutoff scale A, hence leading to a more trustworthy 
perturbative loop expansion of the potential. The first option is to take the SM two-loop 
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result and apply it to each of the quartic couplings at the cutoff scale individually: 



A (A) < Afp/4 
as (A) < Afp/4 
h (A) < Afp/4 



(16) 



where Afp/4 ~ 18 is taken to be a fixed number (in the normalization of Ref. [21], Afp — 12). 
This is the less restrictive of the two perturbative criteria we impose. The second, more 
restrictive choice echoes the constraints imposed by the authors of Ref. j23| (adapted to our 
normalizations) : 

A (A) < An 

a2 (A) < 27r (17) 
h (A) < 2n/3 . 

From a numerical study of the RG-improved one-loop effective potential, we find that using 
either of these two options allows us to avoid the unstable minima in the potential that arise 
when the triviality criterion is used. Additionally, these perturbativity restrictions indicate 
a range where the upper bounds of a more extensive two-loop analysis may lie. 

With these considerations in mind, we translate the bounds on X{Mz) into bounds on 
Mh by working at the scale /i = f o where we impose the minimization condition in Eq. (ITB]1 . 
The corresponding eigenvalues of the mass-squared matrix are given by 



cff 



h=vo, S=0, ii=vo 



cff 



(18) 
(19) 



h=vo, S=0, n=vo 

The unbroken Z2 symmetry of the potential implies that the off-diagonal term 



m 



cff 



hS 



dhdS 



= 



(20) 



h=vo, 5=0, ti=V() 

SO that the doublet and singlet fields do not mix. 

The physical masses of the Higgs and singlet, M| and M|, are determined by the poles 
of the propagators. In the MS scheme, they are related to the running masses mKfi) and 
m|(/i) by solving 



= Ml- (mlif^) + {p' = Ml /x)) 
= M| - (m|(/i) + £5 (p2 = Ml /i) 



(21) 
(22) 

where E/^ 5 denote the MS renormalized self energy functions. The eigenvalues and m| 
appearing in Eqs. f lTSj) and f lT9|) . however, include contributions to the self-energy functions 
evaluated at = 0, since the one-loop effective potential contains the sum of all one-particle 
irreducible one-loop graphs with zero external momentum. Consequently, when relating 
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m\ g to the pole masses M| ^ we must subtract out the = self energy contributions as 



discussed in Ref. 13 Ol: 



= Mi 



m, 



where 



AE, = % (p2 = Mf) - % = 0) 



(23) 
(24) 

(25) 



The self-energy correction to the Higgs AS/j is calculated in the Standard Model in 
the appendix of Ref. [sO] (the authors use the notation AH), but there is an additional 
contribution to the self-energy from the scalar singlet which must be calculated. The singlet 
contributes only one diagram to the Higgs self-energy, and the result is 



AS 



{s) 



Svr 



2^2 {vofvl 



(26) 



where the function Z is defined in Ref. [30 

y 



as 



Z{x)=2+ I dylog 







X 



2y/\l - 4x| arctan - 4a;|j 

V|l -4x|log 



1 + 



4x| ) / (1 - V|l -4x| 



for a; > - 
4 

forx < - 
4 



On the other hand, there is only one diagram which contributes to AS5. This diagram gives 
the result 



AS. 



1 . ,2 2 f\ . Ml{l~x) + Mlx^ 
ci2{vQ) dxlog- 

Jo 



47r2 
1 

4^ 



Ml (1 - a;) + Mix 



(27) 



02 {vofvl 



Ml 



Ml 



Ml Ml\^ 

— loe — ^ - 1 

Ml -Ml 2Ml) ^Ml . 



where the new function y is defined as 
1 

y{x) 



X 



V4x - 1 



arccot (V4x - 1) + arctan 



2a;- 1 
V4x - 1 



X 



Vl — 4x arctanh (-\/l — 4a;) 



for a; > 
forx < 



(28) 



The physical masses can now be calculated from Eqs. ( 123|1 and ( l2l|l . 



4. RESULTS & DISCUSSION 



From the expression for /3a in Eq. f ICSI) . we observe that the presence of the quartic 
WHS^ operator generates a positive contribution to the running of A. One would thus 
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expect that the onset of a turnover in Ves along the h direction to occur at larger scales 
than in the Standard Model, while the occurrence of the perturbativity bound would occur 
at relatively lower scales. As a result, both the vacuum stability and triviahty bounds 
on Mh in the Z2XSM should be lowered for a given A, relative to those in the SM. This 



effect was demonstrated in Refs. [23|, l2J] for A Mpi. In what follows, we verify this 
expectation explicitly for cutoffs from 1 TeV to 10^^ GeV ~ Mpi and for suitable choices of 
the boundary conditions on 62 and 64 leading to minimal impact of these parameters on the 
vacuum structure of the theory. 

In general, however, the singlet quadratic and quartic parameters can play a non-trivial 
role. To see why, consider Eq. ([5]). Depending on the value chosen for a2, a negative value 
for 62 at the input scale may be required to obtain a desired value of m|. A negative value 
for 62 (a^o) Kiay lead to an additional minimum for 5* 7^ that will be deeper than the 
electroweak minimum unless 64 is sufficiently large. These considerations are particularly 
important for a scalar dark matter candidate that is relatively light. The coupling 02 controls 
the annihilation cross section, and it must be sufficiently large to prevent a relic density that 
is too large. On the other hand, the corresponding value of the light mass may require a 
negative 62- In addition, larger values of 64 at the input scale can cause both 64 and 02 to 
reach their perturbativity bounds sooner, modifying the perturbativity upper bounds based 
solely on the running of A. Below, we analyze these effects in detail^. 

We first illustrate the impact of the new interactions associated with the singlet scalar 
on the stability of the electroweak minimum and the occurrence of the perturbative limits 
in the Higgs quartic self-coupling. In Fig. [T] we show the range of allowed Higgs masses as 
a function of its coupling to the singlet, 02, while requiring stability of the potential up to 
10^^, 10^, 10^, or 10'^ GeV. In Fig. |2]we plot against the cutoff scale for representative 
values of 02. For Figs. [1] and [21 we have taken the somewhat artificial choice b2{vo) = 0, 
64 (M^) = 10"'^ to isolate the dependence of the Higgs mass bounds on the 02 parameter. 
In Fig. [21 the curves for successively larger values of 02 illustrate the reduction in both the 
stability lower bound and perturbativity upper bound expected from the behavior of /3a- 
For sufficiently large 02, the impact of including this parameter in the one- loop /3 function 
for the Higgs quartic coupling leads to a maximum A (see the orange and purple curves) 
- a feature that does not arise in the SM or for very small values of 02 [e.g., the green 
curve) In short, for a sufficiently strong H'^HS^ interaction, it is not possible to choose a 
value of Mh that maintains both electroweak vacuum stability and perturbative behavior to 
arbitrarily large scales. 

The solid-colored and striped regions in Fig. [1] indicate the choice of perturbativity re- 



^ Note that we only consider positive values of the coupling 02. In principle, a negative coupling is allowed, 
in which case stability of the electroweak minimum imposes a tree level constraint on a2 ■ When including 
quantum corrections and one-loop running, Eq. (jC6|) indicates that a negative value at the input scale 
would cause 02 to become even more negative very quickly with increasing scale (due largely to the top 
quark contribution). This would cause the effective potential to be unbounded below in some direction of 
field space between the h and S axes. Requiring stability of the potential up to a cutoff scale A would then 
impose a lower bound on 02- A negative 02 coupling can contribute to a strong first-order electroweak 
phase transition [llj] , but does not affect the dark matter phenomenology because it enters the equations 
for the annihilation and direct detection cross sections as a\. 

^ In the text, we refer to features of our figures according to their coloring. In the captions of the figures, 
we provide both color and grayscale schemes. A full color version of this work is available online. 
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FIG. 1: The change in the Higgs mass upper and lower bounds as 02 (Mz) is increased while 
requiring stability of the potential and perturbativity of the quartic coupling constants up to cutoff 
scales A = (lO^^, 10^, 10^, 10^) GeV (the purple [dark gray], blue [light gray], green [medium 
gray], and orange [light gray] regions, respectively). Here, 64 (M^) = 0.001 and 62(^^0) = 0. 
Perturbativity of the potential using the constraints of Eq. ()17p is indicated by the solid-colored 
regions; the combined solid-colored and striped regions indicate perturbativity of the potential 
using Eq. ([16]). 



strictions: the solid-colored regions are perturbative up to the indicated cutoff scale when 
the more restrictive choice of Eq. f[T7|) is applied, while the combined solid and striped re- 
gions of each color use the less restrictive choice of Eq. f lT^ to maintain perturbativity to 
the desired A. The width of the upper bands in Fig. |2] show the variation in the Higgs 
mass upper bounds when varying the perturbativity restrictions for the quartic couplings 
at the cutoff scale between Eqs. (fT7|) and (fT6|) . as discussed in the previous section. We 
have not shown any theoretical ranges for these bounds due to other considerations, e.g., 
varying the top quark mass over its experimental range. Compared to the expected impact 
of including two-loop corrections and the ambiguities in choice of perturbativity bounds, the 
Mt uncertainty is subleading. 

In Figs. OandlU we show analogous results to Figs. [T] and [2] where 64 (Mz) is now taken 
to be 0.4 (62 (vq) is still zero). The lower bounds from stability of the potential in the h 
direction are largely unchanged. On the other hand, the upper bounds from perturbativity 
are reduced. The larger value of 64 can cause the couplings to reach their perturbative limits 
below the desired cutoff scale, and so either A (Mz) must be reduced or it is impossible to 
maintain perturbativity of the couplings up to the desired cutoff. In particular, the green 
curve of Fig. IHnow demonstrates maximum scale behavior similar to the orange and purple 
curves, so no value of and 02 will give a stable potential with perturbative couplings up 
to Mpi. 
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FIG. 2: One-loop vacuum stability and perturbativity bounds on the Standard Model Higgs mass 
Mfi as a function of the cutoff A for a set of values of 02 (M^). We show a2 {Mz) = (0.005, 0.4, 0.8) 
using green [medium gray], orange [light gray], and purple [dark gray], respectively. The upper 
bounds show the variation between the less restrictive constraints of Eq. ()16p and the more restric- 
tive constraints of Eq. (fT7|) . 

The impact of including non-zero values for both the singlet quadratic and quartic cou- 
plings can also be significant, particularly when considering the relic density of singlet dark 
matter'^. The 02 and Ms parameters control the SS annihilation cross section and, thus, the 
singlet relic density, assuming thermal production j2|-0]- As discussed in Section [3l the region 
of relatively light Ms and larger 02 can force the parameter 62 ("^^o) to be negative, leading 
to a secondary extremum with 5* 7^ 0. Thus, for some choices of the Z2XSM couplings and 
masses, the extrema with S* 7^ is a minimum of lower energy than the electroweak/dark 
matter (EW/DM) minimum corresponding to [h = vq, S = 0); for other choices of the 
parameters, the extrema with 5* 7^ either do not exist or have greater energy than the 
desired EW/DM minimum. 

To illustrate, we show in Fig. [5] a plot of the coupling 02 vs. Ms- We have fixed the Higgs 
mass, Mh = 200 GeV, and the singlet quartic self-coupling at the input scale, 64 (Mz) = 0.4. 
A cutoff scale of A = 1 TeV is used. The red curve in Fig. O marks the boundary between 
the two regions discussed in the previous paragraph: for values of 02 (Mz) and Ms located 
in the lower unshaded region, the effective potential has a stable EW/DM minimum, while 
02 and Ms values in the shaded upper region are excluded due to the presence of unviable 
S" 7^ minima along the S direction of the potential below the cutoff scale. This boundary 
is the result of our theoretical vacuum stability analysis. 



For simplicity, we henceforth use only the perturbativity constraints in Eq. ()16l) . 
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FIG. 3: The analogous plot to Fig. [T]for the case 64 (Mz) = 0.4. Here there is no region of Higgs 
masses and S-H couplings 02 that allow for perturbativity of the potential up to 10^^ GeV. The 
color and grayscale schemes are the same as Fig. [TJ 
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FIG. 4: The analogous plot to Fig. [2] for the case 64 {Mz) = 0.4, demonstrating that stability 
and perturbativity of the potential cannot be maintained to arbitrarily large {i.e., O {Mpi)) scales 
when 64 is increased. The color and grayscale schemes are the same as Fig. [2j 
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Ms (GeV) 

FIG. 5: Values of 02 (Mz) and Ms consistent with WMAP relic density measurements (blue [dark 
gray] band) and vacuum stability/perturbativity (unshaded regions) for fixed Mf^ = 200 GeV and 
64 (Mz) = 0.4. The unshaded region below the red [medium gray] curve has stable {h = vq, S = 0) 
minima up to a A = 1 TeV cutoff scale. The shaded region above the red curve is excluded because 
of the occurrence of deeper minima along the S direction of the effective potential. Additionally, 
we indicate the region excluded by CDMS-II (brown [dark gray] arrow) and the region of projected 
sensitivity for SuperCDMS (green [light gray] arrow). 



We also show in Fig. various features relevant for dark matter. Given the model 
parameters M^, Ms, and 02, we calculate the annihilation cross section (Xann and relic density 
Qsh'^ for singlet cold dark matter (CDM), following the discussions in Refs. [sl, The blue 
band in Fig. [5] shows the regions required to obtain a relic density consistent with the WMAP 
range^°, assuming that the singlet saturates the CDM relic density. We also calculate the 
cross section for direct detection of singlet dark matter by elastic scattering off a nucleon, 
add, according to Refs. [gI, |18[; this cross section is then scaled by the fraction of the relic 
density constituted by the singlet, 

Qsh"^ ^sh^ 

O'dd — O-dd — CTdd TTTTiT < CTdd , 

ildmh^ 0.105 

as in Ref. to account for reduced flux in detector experiments in scenarios where the 
singlet does not saturate the CDM relic density. The computed ddd is then compared to 



^ We have also made use of the Higgs decay width program HDEC AY of Ref. [35j . 

° These bands in Figs. El El and|9]are reproduced from Ref. [6j] using the 90% C.L. range for the relic density 
taken from the Particle Data Group 0.092 < fid^/i^ < q.hs. 
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FIG. 6: Values of 02 (Mz) and Ms consistent with WMAP relic density measurements (blue 
[dark gray] band) and vacuum stability/perturbativity (unshaded regions). The rows correspond 
to Mh = 120 GeV (top), 200 GeV (middle), and 350 GeV (bottom). The columns correspond to 
64 (Mz) = 0.001 (left), 0.4 (middle), and 0.8 (right). Shaded regions above the red [medium gray] 
curves are excluded because of deeper minima along the S direction. The unshaded regions have 
stable EW/DM minima for a cutoff scale of A = ITeV. Additionally, we indicate the regions 
excluded by CDMS-II (brown [dark gray] arrows) and the regions of projected sensitivity for 
SuperCDMS (green [light gray] arrows). 



the experimental bounds^^ of XENONIO [38| and CDMS-II |39| as well as the projected 



sensitivity of SuperCDMS [17|. We show in Fig. |5]the regions of (02, Ms) parameter space 
excluded by CDMS-II (brown arrow) and the regions to which SuperCDMS will be sensitive 
(green arrow). We do not show the exclusions from XENONIO because they are superseded 
by the more restrictive CDMS-II limits. 

We plot the vacuum stability regions for a set of representative values of Mh and 64 (Mz) 
in Fig. O Each row of plots corresponds to a single Higgs mass (top row, Mh = 120 GeV; 
middle row, Mh = 200 GeV; bottom row, Mh = 350 GeV) while each column corresponds to 
a single value of the singlet self-coupling at the input scale (left column, 64 (Mz) = 0.001; 
middle column, 64 (Mz) = 0.4; right column, 64 (Mz) = 0.8). A cutoff scale of A = 1 TeV 
is used. The red curves of Fig. |6] mark the boundary between the regions with stable 



See Ref. [37[ for a compilation of dark matter detection limits from many different sources. We thank 
M. McCaskey for providing numerical dark matter detection limits as a function of the dark matter mass. 
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FIG. 7: The minimum singlet mass compatible with both the WMAP relic density constraints 
(assuming the singlet saturates the relic density) and vacuum stability/perturbativity as a function 
of the Higgs mass, for three representative values of the singlet quartic self-coupling 64 {Mz)- The 
cutoff scale is taken to be 1 TeV. 

EW/DM minima (unshaded lower areas) and regions which are excluded due to the presence 
of unviable S* 7^ minima (shaded upper areas). The blue bands, as in Fig. |5|, outline the 
regions consistent with the observed relic density, while the brown and green arrows indicate 
the regions of CDMS-II exclusion and SuperCDMS sensitivity, respectively. 

Fig. E] demonstrates that 64 (M^) must be sufficiently large to have a stable elec- 
troweak/dark matter minimum for a given Ms, 02, and M^. In other words, increasing 
64 (Mz) from the left plots to the right plots of Fig. [6] permits larger values of 02 for a given 
Ms and Mh. There are two implications of this. First, for sufficiently small 64, one may 
encounter a minimum value of Ms - given by the intersection of a given line with the blue 
band - for which the {h = vo, S = 0) vacuum is stable and for which S saturates the CDM 
relic density. Increasing the Higgs mass for a fixed 64 (Mz) also increases the minimum 
singlet mass allowed by vacuum stability, as evidenced by moving top to bottom down a 
column of plots in Fig. El We show this effect in more detail in Fig. [7] where we plot the 
minimum M5 as a function of the Higgs mass for the same representative values of 64 (Mz). 

However, we see from Fig. that except in cases of larger M^ and very small 64 (Mz) the 
experimental limits on the singlet mass are more restrictive than the combined theoretical 
vacuum stability and relic density limits shown in Fig. [71 For a relatively low cutoff scale, 
i.e. A = 1 TeV, it is unlikely that discovery of both a Higgs and dark matter scalar, and 
measurement of their masses, would divulge information on the S self-coupling by way of 
Fig. [71 in which the 5* is assumed to saturate the relic density. 

The second implication is the existence of vacuum stability bounds on the singlet relic 
density as a function of the quartic coupling 64 for scenarios in which the singlet is no 
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FIG. 8: The minimum singlet relic density, as a fraction of the total relic density, vs. 64 (Mz) for 
Ms = 65 (purple [dark gray] curve), 80 (orange [light gray] curve), and 95 GeV (green [medium 
gray] curve) and taking Mh = 120 GeV and A = 1 TeV. The regions below the respective curves are 
disallowed because they correspond to larger 02 values which are excluded when requiring vacuum 
stability up to A = 1 TeV. 



longer required to saturate the relic density. Fig. |6] shows that there is a maximum 02 
- corresponding to the red curves - allowed by vacuum stability when keeping the other 
parameters {Ms, Mh, (Mz) , A) fixed. The singlet relic density is inversely proportional 
to the annihilation cross section which itself is proportional to thus, the maximum 

value of 02 allowed by vacuum stability gives a minimum relic density for the singlet. We 
plot in Fig. [H] this minimum S relic density, as a fraction of the total observed relic density 
{VLdmh'^ = 0.105 [sgI), vs. 64 {Mz). We use a cutoff scale of 1 TeV, the Higgs mass is fixed 
at 120 GeV, and we choose three fixed values for the singlet mass, Ms = 65 (purple curve), 
80 (orange curve), and 95 GeV (green curve), all of which fall in the region not excluded 
by CDMS-11 but accessible by SuperCDMS (see the top row of Fig. [6]). We observe that a 
scalar singlet with a mass in the vicinity of 65 — 80 GeV and constituting 2% of the relic 
density would require a quartic self-coupling 64 {Mz) > 0.2 to be consistent with vacuum 
stability. Measurement of Ms and the singlet's coupling with the Higgs (or determination 
of its relic density by other means) could therefore provide enough information to constrain 
the singlet's self-interaction features from theoretical considerations of vacuum stability if S 
is only a minor constituent of the total dark matter relic density. 

Finally, we comment on the effect of increasing the cutoff scale. In the left column of Fig. [9] 
we show again the allowed regions of (02, Ms) space with 64 (Mz) = 0.4 (the middle column 
of Fig. |6]) for a 1 TeV cutoff. In the right column of Fig. [9l we show the same parameter 
spaces but now we require stability and perturbativity of the potential up to A = 10^ GeV. 
For Mh = 120 GeV (top right), the (02, Ms) region below 02 — 0.25 is disallowed due to the 
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FIG. 9: The column on the left is the same as the middle column of Fig. [6]- = 120, 200, 350 GeV 
top to bottom and 64 [Mz) = 0.4 with A = 1 TeV; arrows indicate the regions accessible to CDMS II 
(brown [dark gray] ) and SuperCDMS (green [light gray] ) . The right column shows how the allowed 
(unshaded) regions are further constrained by increasing the cutoff scale to A = 10^ GeV. For 
Mh = 120 GeV, intersection with the vacuum stability bounds due to large h minima removes the 
allowed (02, Ms) region below 02 — 0.25. For = 200 GeV, intersection with the perturbativity 
bounds removes the allowed regions above a2 — 0.35. For M/^ = 350 GeV, perturbativity limits 
eliminate the entire region of (02, M^) space. These effects are consistent with Fig. [3l 

occurrence of large h minima. For Mh = 200 GeV (middle right), the perturbativity bounds 
remove the allowed regions above 02 — 0.35, while the perturbativity bounds prohibit all of 
the (02, Ms) parameter space shown for Mh = 350 GeV (bottom right). These restrictions 
are on 02 only since they are independent of the singlet mass for the cases shown; they are 
consistent with Fig. [3] from which we can determine the values of 02 that allow stability and 
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FIG. 10: The dashed hnes show the minimum fractional singlet relic densities from Fig. [8] vs. 
64 (M^) when requiring vacuum stability to ITeV. Increasing the cutoff scale to A = 10^ GeV 
severely restricts the allowed regions of 02, i.e., the singlet relic density, and 64: only the small 
colored regions are allowed due to vacuum stability and perturbativity restrictions. Mh = 120 GeV 
and Ms = 65 (purple [dark gray]), 80 (orange [light gray]), and 95 GeV (green [medium gray]) as 
in Fig. E 



perturbativity of the potential up to a 10^ GeV cutoff once the Higgs mass is fixed^^. 

More restrictive bounds on the relic density of the singlet can be obtained by requiring 
stability and perturbativity of the potential to higher scales. For example, the existence 
of a minimum 02 coupling in the top right plot of Fig. [9] implies a maximum singlet relic 
density less than the total relic density. We show in Fig. [TU] another plot of the fraction 
Qs/^dm VS. 64 (Mz) for Mh = 120 GeV and Ms = (65, 80, 95) GeV (purple, orange, and 
green, respectively), but now we require stability and perturbativity of the potential up to 
A = 10^ GeV. Whereas in Fig. |8]the regions above the curves (all the way up to Qs/^dm = 1) 
were allowed for A = 1 TeV, increasing the cutoff scale to A = 10^ GeV greatly diminishes 
the allowed regions to the small colored areas of Fig. [10] (the 1 TeV restrictions are shown as 
dashed curves in Fig. [10] for comparison). Assuming a Higgs mass of 120 GeV, we observe 
that imposing stability and perturbativity of the effective potential up to a 10^ GeV cutoff 
scale restricts the singlet to scenarios in which it constitutes 2% or less of the CDM relic 
density and has a moderate (0.2 to 0.5 in our normalization) self-interaction coupling. 



Recall that the less restrictive choice of perturbativity constraints corresponds to the combined solid- 
colored and striped regions of Fig. |3j 
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5. CONCLUSIONS 

We have performed a study of the Z2XSM using requirements of vacuum stabihty and 
perturbativity of the one-loop effective potential. The extra scalar degree of freedom intro- 
duced in this model can provide a suitable cold dark matter candidate. In this scenario, 
the relevant parameters, in addition to the Higgs mass, are the singlet mass M5, its 02 
coupling to the Higgs by the WHS^ term of the potential, and the quartic self-coupling 
64. Ms and 02 govern the annihilation cross section of the singlet, its relic density, and its 
cross section for scattering off a nucleon. By imposing the requirements of vacuum stability 
and perturbativity of the potential, we are able to constrain the parameter space of the 
Z2XSM. If the singlet saturates the relic density, stability of the {h = vq, S = 0) minimum 
of the potential up to the 1 TeV scale necessitates a lower bound on M5 that varies with 64. 
Measurement of Ms could therefore in principle constrain the singlet self-coupling; however, 
current experimental limits already exclude a wide range of singlet masses, thus limiting our 
ability to constrain the self-interactions. If the singlet undersaturates the relic density, we 
can again place limits on 64 from vacuum stability; however, for values of Ms not already 
excluded by experiment, such limits are relevant when using a 1 TeV cutoff for our theory 
only if the singlet constitutes O (1%) of the total relic density. Increasing the cutoff scale 
places much stronger bounds on the singlet relic density and its self-coupling. For instance, a 
scalar singlet discovered by SuperCDMS would compose at most 2% of the relic density and 
have a moderate self-coupling in the range 0.2 — 0.5 if we require stability and perturbativity 
of the effective potential up to 10^ GeV. 

These analyses of the Z2XSM parameter space just described depend on the as-yet undis- 
covered Higgs mass; we have also studied the reverse effect in which the Z2XSM parameters 
affect bounds on the Higgs mass arising from vacuum stability and perturbativity. Our nu- 
merical studies found that triviality, i.e., avoidance of Landau poles below the cutoff scale, 
is an insufficient criterion for determining upper bounds on the Higgs mass, forcing us to 
require perturbativity of the running couplings instead. Our analysis showed that increas- 
ing the 02 coupling decreases both vacuum stability and perturbativity bounds on the Higgs 
mass, while increasing 64 decreases only the upper bounds from perturbativity. Additionally, 
we have shown that requiring vacuum stability and perturbativity of the potential for larger 
02 and 64 can force the appearance of new physics at a scale below the Planck scale. 

In the case of the Z2XSM, theoretical requirements of vacuum stability and perturbativity 
have interesting implications for the relationships between the singlet scalar self-interaction, 
the corresponding fraction of the dark matter relic density, and the scale at which one expects 
additional degrees of freedom to become relevant. These considerations can also affect the 
theoretical interpretation of searches for the SM Higgs boson and for scalar dark matter in 
direct detection experiments. It is natural to ask if there exist similar implications for other 
non-singlet scalar dark matter scenarios, such as those in which additional scalar degrees 
of freedom have non-trivial gauge quantum numbers (e.g., the triplet model discussed in 
Ref. [4^). Also of interest are scenarios in which extra scalars can contribute to a strong 
first-order electroweak phase transition. Analyses for these situations are the subjects of 
ongoing work. 



22 



Acknowledgements 

This work was supported in part by Department of Energy contract DE-FG02-08ER41531 
and the Wisconsin Alumni Research Foundation. We thank P. Langacker, D. J. H. Chung, 
and D. Hooper for useful discussions and M. McCaskey for providing the numerical dark 
matter detection limits as a function of the dark matter mass. We thank M. Sher for a 
careful reading of this manuscript and useful comments. MJRM also thanks NORDITA, 
the Aspen Center for Physics, and the TRIUMF Theory Group where part of this work was 
completed. 



Appendix A: Notation & Normalization of the Potential 

Here we compare our normalization choices for the parameters of the potential, Eq. ([T]), 
with other choices found in the literature. In the convention commonly found in the literature 



{e.g. Refs. l2l|, |23|, l2J|) where ml = 2Xvq, the coupling A is a factor of six smaller than 
our coupling A. Our normalization choice of m\ = Xvq/3 follows that of Ref. [2^ which 
was our source for the SM RGEs; this choice leads to a Higgs quartic self-coupling term 
in the potential preceded by ^, thus mimicking conventional 0^ theory. Additionally, our 
Higgs-singlet couplings Oi and 02 are a factor of two larger than those found in Refs. pjj], but 
smaller by a factor of two compared to Ref. (in which the authors use the notation X^s 
rather than 02 and do not include an ai type term because of their choice of a Z2 symmetry 
from the outset). 

Appendix B: One-loop Potential Particle Content 

The numerical factors and field- dependent masses in Eq. ([6]) are as follows for the indi- 
cated particles: 
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m\ = see Eqs. §^ and ([7]) 



Here, stands for the would-be Goldstone bosons and ± stands for the diagonalized 
Higgs/singlet contributions. We have not included contributions to the one-loop potential 
from fermions other than the top quark due to the smallness of their Yukawa couplings. 
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Appendix C: One-loop Beta Functions 

We give here the /3 and 7 functions for all parameters of the Z2XSM. All couplings depend 
on the 't Hooft scale n which we have suppressed for clarity. 
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